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■ Abstract 

We study the resummation of large logarithmic perturbative corrections to the partonic cross 

•rH , 

sections relevant for the process pp — ► W^X at the BNL Relativistic Heavy Ion Collider (RHIC). 
At RHIC, polarized protons are available, and spin asymmetries for this process will be used for 
precise measurements of the up and down quark and anti-quark distributions in the proton. The 
corrections arise near the threshold for the partonic reaction and are associated with soft-gluon 
emission. We perform the resummation to next-to-leading logarithmic accuracy, for the rapidity- 
differential cross section. We find that resummation leads to relatively moderate effects on the 
cross sections and spin asymmetries. 
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I. INTRODUCTION 



The exploration of the spin structure of the nucleon has by now become a classic topic 
of nuclear and particle physics. Most of our knowledge so far comes from spin-dependent 
deeply-inelastic scattering (DIS) which has produced spectacular results Q, the single most 
important one being that the total quark spin contribution to the nucleon spin is only about 
20 — 25%. This has led to a new generation of experiments that aim at further unraveling 
the spin structure of the nucleon. Among these are experiments at the Relativistic Heavy 
Ion Collider (RHIC) at BNL. Here, a new method is used, namely to study collisions of 

n 

two polarized protons 2]. This allows in particular clean probes of the polarized gluon 
distribution in the proton, whereby one is hoping to learn about the contribution of gluons 
to the proton spin. 

DIS has given us exciting information on the spin-dependent quark and anti-quark 
distributions Aq(x,fi 2 ), Aq(x,fi 2 ), where Aq(x,fi 2 ) = q\x,fi 2 ) — (^(x, u 2 ) with q^(x, « 2 ) 
(q^(x,fi 2 )) denoting the distributions of quarks with positive (negative) helicity and light- 
cone momentum fraction x in a proton with positive helicity, at factorization scale «. In- 
clusive DIS via photon exchange, however, gives access only to the combinations Aq + Aq. 
In order to understand the proton helicity structure in detail, one would like to learn about 
the various quark and anti-quark densities, Am, Am, Ad, Ad, As, As, individually. This is 
for example relevant for comparisons to models of nucleon structure which generally make 
clear qualitative predictions about, for example, the flavor asymmetry Am — Ad in the pro- 
ton sea Q|. These predictions are often related to fundamental concepts such as the Pauli 
principle: since the proton has two valence-M quarks which primarily spin along with the 
proton spin direction, mm pairs in the sea will tend to have the u quark polarized opposite 
to the proton. Hence, if such pairs are in a spin singlet, one expects Am > and, by the 
same reasoning, Ad < 0. Such questions become all the more exciting due to the fact that 
rather large unvolarized asymmetries u — d ^ have been observed in DIS and Drell-Yan 
measurements 0, 0, |f| . 

There are two known ways to distinguish between quark and anti-quark polarizations 
experimentally, and also to achieve at least a partial flavor separation. Semi-inclusive mea- 
surements in DIS are one possibility, explored by SMC |7j and, more recently and with higher 
precision, by Hermes jSj]. Results are becoming available now also from the Compass ex- 
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periment P| , and measurements have been proposed for the Jefferson Laboratory . One 
detects a hadron h in the final state, so that instead of Aq + Aq the polarized DIS cross 
section becomes sensitive to Aq(x) D^(z) + Aq(x) D^(z), for a given quark flavor. Here, 
the D*l(z) are fragmentation functions, with z denoting the fraction of the momentum of 
the struck quark transfered to the hadron. The dependence on the details of the fragmen- 
tation process limits the accuracy of this method. There is much theory activity currently 
on SIDIS, focusing on next-to-leading order (NLO) corrections, target fragmentation, and 
higher twist contributions 111 ]. 

The other method to learn about the Aq and Aq will be used at RHIC. Here one considers 
spin asymmetries in the production of W ± bosons [12J. Because of the violation of parity in 
the coupling of W bosons to quarks and anti-quarks, smg/e-longitudinal spin asymmetries, 
defined as 



da ++ + dcr + ~ - d(J~ + - dor" dA 



a 



Al da++ + da+~ + da~+ + do— do- ' ^ 
can be non- vanishing. Here the a denote cross sections for scattering of protons with definite 

helicities as indicated by the superscripts; as one can see, we are summing over the helicities 

of the second proton, which leads to the single-spin process pp — > W X. Since the W mass 

sets a large momentum scale in the process, the cross sections factorize into convolutions 

of parton distribution functions and partonic hard-scattering cross sections. The latter 

are amenable to QCD perturbation theory, the lowest-order reaction being the "Drell-Yan" 



process qq' — > W . For produced W + , the dominant contri 



and the spin asymmetry is approximately given by 



Bill [3 



3utions come from ud — > W + , 



w+ = Au(x° 1} M^)d(x° 2} Mir) - Ad(x° 1} M^)u(x° 2} M^) 

L ~ u{xlM^)d{xlM" w ) + d{xlM^)u{xlM^) ' l ' J 

Here, as indicated, the parton distributions will be probed at a scale ~ M w , with M w the 

W mass, and the momentum fractions x\ and x\ are related to the rapidity 77 of the W + by 

< 2 = ^p', (1-3) 

where \fS is the pp center-of-mass energy which at RHIC is 200 or 500 GeV. Note that rj 
is counted positive in the forward region of the polarized proton. From Eq. (J1.2j) it follows 
that at large 77, where x\ ~ 1 and x® <C 1, the asymmetry will be dominated by the 
valence distribution probed at and hence give direct access to Au(x1, M^) / u(x1, M^r) . 
Likewise, for large negative rj, is given by — Ad(xi, Myy)/d(xi, M^). For negatively 



charged W bosons, one finds correspondingly sensitivity to Ad(x®, Myy)/d(x1, Myy) and 
— Au(x1, Myy)/u(x\, M$y) at large positive and negative W rapidities, respectively. This is 
the key idea behind the planned measurements of Au, Au, Ad, Ad at RHIC. 

In practice, a significantly more involved strategy needs to be used in order to really relate 
the single-longitudinal spin asymmetries to the polarized quark and anti-quark densities. 
Partly, this is an experimental issue: the detectors at RHIC are not hermetic, which means 
that missing-momentum techniques for the charged-lepton neutrino (lv{) final states cannot 
be straightforwardly used to detect the W and reconstruct its rapidity. There are, however, 
workarounds to this problem. One is to assume that the W is dominantly produced with 
near-zero transverse momentum, in which case one can relate the measured rapidity of the 
charged lepton to that of the W , up to an irreducible sign ambiguity \lL\- Ultimately, it 
may then become more expedient not even to consider the W rapidity, but to formulate all 
observables directly in terms of the charged-lepton rapidity. This approach has been studied 
in detail in . It was found that despite the fact that the interpretation of the spin 

asymmetries in terms of the An, Au, Ad, Ad becomes more involved, there still is excellent 
sensitivity to them. Contributions by Z bosons need to be taken into account as well. Very 



recently, it has also been proposed 17[ to use hadronic W decays in the measurements of 
Al. These have the advantage that one just needs to look for events with two hadronic 
jets. A potential drawback is that, while the QCD two-jet background cancels in the parity- 
violating numerator of Al, it does contribute to the denominator and therefore will likely 
reduce the size of the spin asymmetries. 

There are also theoretical issues that modify the picture related to Eq. (jl.2j) described 
above. There are for example Cabibbo-suppressed contributions. It is relatively straightfor- 
ward to take these into account, even though one needs to keep in mind that they involve the 
polarized and unpolarized strange quark distributions. More importantly, there are higher- 
order QCD corrections to the leading order (LO) process qq' — > W . At next-to- leading 
order, one has the partonic reactions qq' — > W ± g and qg — > W ± q' . For the unpolarized and 
the single-longitudinal spin cases, the cross sections for these have been given and studied 
in Ref, Q Q and 3 fl fl H H H Q, respective,, In the present paper, we 



will improve the theoretical framework by going beyond NLO and performing a resumma- 
tion of certain logarithmically enhanced terms in the partonic cross sections to all orders in 
perturbation theory. 
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The corrections we are considering arise near "partonic threshold" , when the initial par- 
tons have just enough energy to produce the W boson. Here the phase space available for 
real-gluon radiation vanishes, while virtual corrections are fully allowed. The cancellation of 
infrared singularities between the real and virtual diagrams then results in large logarithmic 
"Sudakov" corrections to the LO qq' cross section. For example, for the cross section inte- 
grated over all rapidities of the W, the most important logarithms (the "leading logarithms" 
(LL)) take the form aj[ln 2fc_1 (l — z)/(l — z)] + at the kth order of perturbation theory, where 
z = M^y/s with s the partonic center-of-mass energy, and a s is the strong coupling. The 
"-(-"-distribution is defined in the usual way and regularizes the infrared behavior z — > 1. 
Subleading logarithms are down by one or more powers of the logarithm and are referred to 
as "next-to- leading logarithms" (NLL), and so forth. Because of the interplay of the partonic 
cross sections with the steeply falling parton distributions, the threshold regime can make 
a substantial contribution to the cross section even if the hadronic process is relatively far 
from threshold, that is, even if M^/S <C 1. If M^/S ~ 1, the threshold region will com- 
pletely dominate the cross section. This is expected to be the case for W production at the 
RHIC energies, in particular at = 200 GeV. Sufficiently close to partonic threshold, the 
perturbative series will be only useful if the large logarithmic terms are taken into account 
to all orders in a s . This is achieved by threshold resummation. 

Drell-Yan type processes have been the first for which threshold resummation was derived. 
The seminal work in 27, 0] has been the starti ng p oint for the developments of related 
resummations for many hard processes in QCD j29(. The resummation for Drell-Yan is 
now completely known to next-to-next-to-leading logarithmic (NNLL) accuracy 2^. In 
this paper we will, however, perform the resummation only at NLL level. Since the strong 
coupling at scale Mw is relatively small, we expect this to be completely sufficient for a good 
theoretical description. In addition, away from the threshold region, one needs to match the 
resummed calculation to the available fixed-order one. In our case of single-spin production 
of W's, this is NLO. A consistent matching of a NNLL resummation would require matching 
to NNLO, which is not yet available. 

As we discussed above, of particular interest at RHIC for determining the spin- dependent 
quark and anti-quark densities are the rapidity distributions of the W bosons. In order to 
present phenomenology relevant for RHIC, we will therefore perform the resummation for 
the rapidity-dependent cross sections. Note that we will concentrate only on the case of the 
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distributions in W^-boson rapidity; from the earlier discussion it follows that for future studies 
it would be even more desirable to consider the distributions in charged-lepton rapidities. To 
treat the resummation at fixed W rapidity, we will employ the technique developed in [31]. 
This entails to use Mellin moments in r = fS of the cross section, as is customary in 
threshold resummation, but also a Fourier transform in rapidity. In [3l|, this method was 
applied to the prompt-photon cross section. For the case of the Drell-Yan cross section 
it simplifies considerably. We note that techniques for the threshold resummation of the 
Drell-Yan cross section at fixed rapidity were also briefly discussed in j^. We verify the 
main result of that paper, and for the first time, present phenomenological results for this 
case. 

Before turning to the main part of the paper, we mention that important work on another 
type of resummation of the W production cross sections at RHIC has been performed in the 



15 



16, 2S 



literature, namely for the transverse-momentum distribution (p?) of the W's 
At lowest order, the W is produced with vanishing transverse momentum. Gluon radiation 
generates a recoil transverse momentum. By a similar reasoning as above, when tends 
to zero, large logarithmic corrections develop in the pt spectrum of the Ws. These can be 



15. 



16 



For 



resummed as well, which was studied for the case of RHIC in great detail in 
the total or the rapidity-differential W cross section we are interested in, pt is integrated, 
and the associated large logarithms turn partly into threshold logarithms and partly into 
nonlogarithmic terms. Therefore, the threshold logarithms we are considering here become 
the main source of large corrections to the cross section. We note that ultimately it would be 
desirable to perform a "joint" resummation of the px and threshold logarithms, as developed 
in Ref. Q. 

The remainder of this paper is organized as follows. In Section [H] we present general 
formulas for the single-spin cross section for W production as a function of rapidity, and 
also discuss the NLO corrections. In Section HTT| we give details of the Mellin and Fourier 
transforms that are useful in achieving threshold resummation of the rapidity-dependent 
cross section. The next two sections provide the formulas for the resummed cross sections. 
In Section fVTl we present our numerical results for RHIC. 
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II. CROSS SECTION FOR W PRODUCTION IN pp COLLISIONS 



The rapidity-differential Drell-Yan cross section for W ± production in singly-polarized 
pp collisions can be written as 
dAcx 



dr] 



r 1 r 1 I M \ 

MJ2 / dxi / dx 2 AV ij [x 1 ,x 2 ,x 1 ,x° 2 , a s (/i 2 ),— ^ Af i (x u fi 2 )f j (x 2 ,iJ, 

M ^2 °ij J d^i J Q dx 2 



h3 

x^Dgg Lti, x 2 , x° v x° 2 , a s (/i 2 ), 



Aq^Xi, n 2 )qj(x 2 , fi 2 ) + Aqi(x u fi 2 )qj(x 2 , \i 2 



AD gq x 1 ,x 2 ,x 1 ,x 2 ,a s (n ), 



M 2 \ 

o _o „ / ,,2\ m w \ \n(„ 



/// 

M 2 s 

^0 „. t..2\ m W 



AG(xi, ^ 



qj(x 2 ,n ) - qj(x 2 ,n 



+D qg x 1 ,x 2 ,x 1 ,x 2 ,a s (n ) 



-Agj(xi,/i 2 ) + Agi(xi,/i 2 ) G(x 2 , \i 



(2.1) 



where the Afi, fj are the polarized and unpolarized parton distributions, respectively, and 
where in the second equation we have written out explicitly the various contributions that 
are possible through NLO. Furthermore, the normalization factor M is given by 

nG F M w V2 



3S 



with the Fermi constant GV, and 0% j £11*6 the coupling factors for W ± bosons, 



\V--\ 2 



(2.2) 



(2.3) 



with Vij the CKM mixing factors for the quark flavors \i ~ My/ denotes the renormal- 
ization/f actor ization scales which we take to be equal. x® and x\ have been defined above 
in Eq. (JEEJ). 

The Dgg, AD gq) D qg in Eq. (|2.1j) are the hard-scattering functions. For the single-spin 
cross section, always one initial parton is unpolarized. Therefore, thanks to helicity con- 
servation in QCD and the V — A structure of the Wqq' vertex, when the polarized parton 
is a quark, one finds that the spin-dependent partonic cross section is the negative of the 
unpolarized one. We have therefore omitted the A's in these cases and only included a A for 
the case of an initial polarized gluon where there is no trivial relation between the polarized 
and the unpolarized cross sections. Each of the hard-scattering functions Dij (or ADij) is a 
perturbative series in the strong coupling a s ([i 2 ): 



M 2 ' 

™0 „ l,,2\ m W 



Uij I 2-2; 2-1) X 2 i CK S (// 



^ij (^lj %2-> ^2^j 



2vr 



-D{j I 3^1, X 2 , %li 3^2 



/If 2 " 



+ 



Only the gg' process has a lowest-order contribution: 



(2.4) 



Dy(xi,x 2 ,x°,xl) = 8(xi - xl) 5{x 2 - xl) , 

ad;;(x 1i3;2 ,x^2) = ^(m,^,^;,^) = o 



(2.5) 



The NLO hard-scattering functions have rather lengthy expressions which have been derived 



in Refs. 



19 



21 



24j . For the reader's convenience we collect them in the Appendix. As can 



be seen from Eqs. ()2.5|) . (|Al|) - (|A5jl . the coefficients contain distributions in (x\ — x®) and 
(x2 — x®)- These are the terms addressed by threshold resummation. More precisely, it turns 
out that only products of two "plus" -distributions, or a product of a "plus" -distribution and 
a delta-function, are leading near threshold. This occurs only for the coefficient D q q. To be 
able to write down the resummed expressions for arbitrary rapidity, we need to take suitable 
integral transforms, to which we shall turn now. 



III. MELLIN AND FOURIER TRANSFORMS, AND THRESHOLD LIMIT 

Consider the production of a W boson through a single generic partonic reaction involving 
initial partons i and j. We define a double transform of the cross section, in terms of a Mellin 
transform in r = Myy / S and a Fourier transform in the rapidity of the W boson: 

Aa(N, M) = / drr N - 1 / ^ cfye^qp (3.1) 
Jo ^~ ln 7? ""q 



= jVf dx 1 x^ +iM/2 Af i (x 1 ) ! dx 2 x 2 "- lM/2 f j (x 2 ) [ dzz 1 "' 1 1^ dfje^AV^. 

JO Jo Jo i-lni 

Here we have suppressed the argument of the partonic cross section as well as the scale 
dependence of the parton distributions. We have introduced parton level variables z = 
r jx\x 2 and r\ — T] — -^ \sy{x\jx<i). Each of the functions AT>ij is of the form {x\x 2 )~ 1 times a 
function of x\jx\ and x 2 /x 2 only. Defining the moments 

f N Q?) = f'dxx^fix,^), (3.2) 
Jo 

of the parton densities, and 

Af>ij(N, M) ee f 1 dzz*' 1 t~ z dfie iMf > (ziz 2 APy) , (3.3) 

Jo J—\n-^= 



we therefore have: 



Aa(N,M)=NAf l N+iM/2 ff- iM/2 AV ij (N,M) . 



(3.4) 



At lowest order 



■>■]_, x 2 , x®, xfj = Xix 2 5(xi — x®) 8(x 2 — x®) = 8 (l — y/ze^ 5 (l — y/ze ^ , (3.5) 



and we find 



dz z 



,ln4- 
N-l I 



In -4. 



\ dfj e iMf > 8(1- ^) «J (l - v^e 



X dzz N ~ x ~ ' 



2 cos (Mln— j (5(1-^) 
(3.6) 

As one can see, a (5(1 - z) emerges, as expected from the well-known expression for the 
rapidity-integrated LO Drell-Yan cross section 18(. Beyond LO, a closed calculation of the 
double moments AT>ij(N, M) is very difficult. For example, at NLO, one would need to use 
the expressions given in Eqs. (|Alj) - (|A5J) and take the moments in A and M. Fortunately, 
a great simplification occurs in the near-threshold limit. Equation ()3.6|) shows that the 
partonic threshold, after taking Fourier moments in rapidity, is reached at z — 1. In Mellin- 
moment space, this corresponds to the large- A limit. Even though the Cosine factor is 
obviously unity in conjunction with the term 8(1 — z), it is generic, as we will now see. 
At NLO, one finds from Eqs. (|Alj) - (|A5J) the following structure near threshold: 



dz z 



N-l 



In -4= 



j df)e tMf > x 1 x 2 D { g 'q > I x 1 ,x 2 ,x { l,x 2 ,a s (/j 



o jo 



111 



.2, 



^dzz N ' 1 I 2 cos (Mln-^=\ C F 



+ 



+ 



2tt 2 



+ 31n- 



w 



8(1 -z) 



+ 0(l-z) 



z + 



In # 



(3.7) 



V 3 

where Cp = 4/3. The factor in square brackets is the large- z limit of the well-known ^| 
0(a s ) QCD correction to the rapidity-integrated Drell-Yan cross section through the qq 
channel. As usual, the "plus" -distributions are defined over the integral from to 1 by 



1 dzf(z) [g(z)} + = dz (f(z) - f(l)) g(z) . 
10 Jo 

As before, the Cosine factor emerges, which is subleading near threshold since 

/ 1 \ (l-z) 2 M 2 



(3.8) 



cos M In 



O KI - zYM 1 



(3.9) 



It is therefore expected to be a good approximation to set this term to unity, and in any 
case consistent with the threshold approximation. On the other hand, the term carries 
information on rapidity through the Fourier variable M. At very large M keeping the 
Cosine term may be more important. In the following, we will ignore the term and just 
point out what changes to our formulas below would occur if one kept it. We will also study 
the numerical relevance of the Cosine term later in the phenomenology section. 

We now take the Mellin moments of the expression in Eq. (|3.7j) . At large N, in the 
near-threshold limit, the moments of the NLO correction become 



dzz N - 1 C F 



o 

= 2C F 
where 



'\n(l-z)\ 4 . Mycr ( 2vr 2 M 2 



w 



, /3 nl ,r rt \ 1 Ml . 27T 2 



2 ln 2 (iV) + ( - - 2 ln(iV) J In — f - 4 + 



+ ( 3 - 10 ) 



N = Ne^ E . (3.11) 

The main result, expressed by Eq. ()3.10|) . is that near threshold the double moments of 
the rapidity-dependent cross section are independent of the Fourier (conjugate to rapidity) 
variable M, up to small corrections that are suppressed near threshold. The iV-dependence 
is identical to that of the rapidity-integrated cross section. This result was also obtained 
in • Near threshold, the dependence on rapidity is then entirely contained in the parton 
distribution functions: as can be seen from Eq. ()3.4j) . their moments are shifted by in- 
dependent terms. As a further source of rapidity dependence, one can keep the Cosine term, 
as discussed above. It is straightforward to keep this term when taking the Mellin moments, 
by writing 

cos (Mm = 1 (z lM ' 2 + z~ iM ' 2 ) . (3.12) 

This will simply result in a sum of two terms of the form (J3.1U|) with their moments shifted 
by ±iM/2. 

We finally recall that contributions from the quark-gluon channels are not equally en- 
hanced near threshold, but are all down by 1/N. They are therefore not subject to resum- 
mation, but will be included in our numerical studies at NLO level, using Eqs. (jA2|) - ()A5l) . 
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IV. RESUMMED CROSS SECTION 



In Mellin-moment space, threshold resummation for the Drell-Yan process results in the 



exponentiation of the soft-gluon corrections. To NLL the resummed formula is given in the 
MS scheme by 0, 3 



dz z 



N-l 



ln 4= / /If 2 " 

^ df)e iMf > Xl x 2 D™ fx 1 ,x 2l xS,4a.(Ai 2 ),^ 



In ^= 



exp 



1 JA - 1 fWtfr dk 



i-C 



A 7 



A(a s (4)) ,(4.1) 



where 



with 



A,(a s 



7T y V 7T 



4 2) + 



67 
18 



7T 



9 



(4.2) 



(4.3) 



Here Nf is the number of flavors and Ca = 3. The coefficient C q collects mostly hard virtual 



corrections. Its exponentiation was shown in 



3J: 



C q a s (n ),ln 



2^ 



7T 



, 2tt 2 3 . 
-4 + — + - In — ^ 



+ (9(a 2 . 



(4.4) 



Eq. (|4.1j) as it stands is ill-defined because of the divergence in the perturbative running 
coupling a s (kji) at kx = Aqcd- The perturbative expansion of the expression shows factorial 
divergence, which in QCD corresponds to a power-like ambiguity of the series. It turns out, 
however, that the factorial divergence appears only at nonleading powers of momentum 
transfer. The large logarithms we are resumming arise in the region 28] £ < 1 — 1/N in the 
integrand in Eq. (|4.1|) . One therefore finds that to NLL they are contained in the simpler 
expression 



M 2 W 



dk 2 



Mir/m kl 



A q (a s {k 2 T ))hi 



Nkj 



+ 2 



m 2 dk 2 



T 



Ml k\ 



A q (a s (k 2 T ))^N 



(4.5) 
of the 



for the second exponent in (|4.1j) . This form is used for "minimal" expansions 
resummed exponent. 

In the exponents, the large logarithms in N occur only as single logarithms, of the form 
a k \n k+1 (N) for the leading terms. Subleading terms are down by one or more powers of 
ln(iV). Knowledge of the coefficients A q 1,2 ^ in Eq. ()4.1|) is enough to resum the full towers 
of LL terms a k ln k+1 (N), and NLL ones a k \n k (N) in the exponent. With the coefficient 
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C q one then gains control of three towers of logarithms in the cross section, a^ln 2fc (iV), 
a k s ln 2k -\N), a k s \n 2k - 2 (N). 

To NLL accuracy, one finds from Eqs. gU),^ HQ 

"1 /-In 4= / M 2 \ 



J dzz N - 1 J ^ dfie iM * Xl x 2 Dlf lx l ,x 2 ,x° l ,x 2 ,a s {^ 



/< 2 



exp{c 9 (a.(A^)+21n^^^ , (4.6) 



where 



A = 6 a s (/i 2 )lniV. (4.7) 



The functions hP~ ,2 > are given by 



h<»(\) = [2A + (1 - 2A) ln(l - 2A)] , (4.8) 



h (2) ( Xi Mk) = -^L- [2A + ln(l - 2A)] + ^ 



2A + ln(l - 2A) + - ln 2 (l - 2A) 
2 



/r J 2n 2 b 2 n 2yr6g 

AW M 2 AWffl (u 2 ) M 2 

\- [2A + ln(l - 2A)] In % - g ^ ' In iV In % , (4.9) 

27TO /i 2 7T /i 2 



where 



60 = -^(11^-2^) , 

61 = ^(17^-5^^-3(7^) . (4.10) 

The function hP~> contains all LL terms in the perturbative series, while is of NLL only. 
When expanded to 0(a s ), Eqs. f)4.6|) - f)4.9|) reproduce the full expression (|3.10|) for the NLO 
correction at large N. We note that the resummed exponent depends on the factorization 
scales in such a way that it will compensate the evolution of the parton distributions. This 
feature is represented by the last term in ([4.9)1 . One therefore expects a decrease in scale 
dependence of the cross section from resummation. The remaining //-dependence in the 
second-to-last term in (14.91) results from the running of the strong coupling constant. 

A s wa S shown in JkQ, it is possible to inwove the above tannla sb g b t l y and t o 
also correctly take into account certain subleading terms in the resummation. To this end, 
we rewrite Eqs. (I4.6|) - (|4.9j) as 

l r lnJ h I M 2 \ 

dzz N ~ l \ * dfje iMf 'x l x 2 D I q f [x l ,x 2 ,x 1 ,x° 2 ,a s (fi 2 ),^} 
(i J ~ h1 ^ \ V J 
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°0 



q n 2 



7T 



4% 

2nb 3 



hl2(1 _ 2A) + B (.»Mi_^) 

7T0 



+ 



' a a {jH 2 ) ^M& | ln(l-2A)' 



7T 



(4.11) 



where M 1 ) = -3CV/2. The last term in Eq. (gUJD is the LL expansion of the term 



M w/* 2 dk\ a s {k 2 T ) 



Since the factor 



-2A<P IniV-SW 



7T 



InN-B^ 



(4.12) 



is the large- A" limit of the moments of the LO splitting 



function APj^, the term in Eq. (|4.12|) may be viewed as the flavor non-singlet part of the 
evolution of the quark and anti-quark distributions between scales u and Mw/N. This 

nn 

suggests to modify the resummation by replacing [38, 39] 



lniV-^W 



C f 



--2Sx(N) + — 

2 U J N(N + 1) 



(4.13) 



in Eq. (|4.11)1 . the term on the right-hand-side being the moments of the full LO non-singlet 
splitting function. With the help of this term, not only the leading large- N pieces of the 
NLO qq' cross section are correctly reproduced, but also the contributions ~ \n(N)/N. We 
note that there are also related h\.(N)/N pieces in the NLO cross section for the qg channel. 
These could be taken into account by extending (J4.13j) to the singlet evolution case, which 
we however refrain from in the present paper. 

V. INVERSE TRANSFORMS AND MATCHING 



The final step in arriving at the resummed rapidity-dependent cross section is to take the 
Mellin and Fourier inverse transforms of Aa(N, M) back to the variables r and r\: 



dAa r 



1 

2^ 



dMe~ iMri 



C+ioo 



dNr~ N Aa res (N, M) 



(5.1) 



drj 2ix J-oo 2iri Jc 

Care has to be taken when choosing the contour in complex N space because of the cut sin 



gularities in the resummed exponent. Adopting the "minimal prescription" of |36[ for the ex- 
ponents, we choose the constant C in (|5.1|) so that all singularities in the integrand are to the 
left of the integration contour, except for the Landau singularity at AT = exp(l/2&oa s (/-t 2 )), 
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which lies to the far right. The contour is then deformed j3] into the half-plane with neg- 
ative real part, which improves convergence while retaining the perturbative expansion. In 
this deformation, we need to avoid the moment-space singularities of the parton densities, 
which are displaced parallel to the imaginary axis by ±M/2, as seen from Eq. (J3.4j) . Thus, 
the intersection at C of the contour with the real axis has to lie far enough to the right that 
the contour does not pass through or below the singularities of the parton densities. The 
technique we use to achieve this is described in detail in Ref. (31 1. 

In order to keep the full information contained in the NLO calculation, we perform 
a "matching" of the NLL resummed cross section to the NLO one. This is achieved by 
subtracting from the resummed expression in Eq. (|5.1|) its 0(a s ) expansion, 



1 roo 1 rC'+ioo 

— / dMe~ lMri — / dNr- N 
2tc J-00 2ni Jc- 



-100 



Ad ies (N, M) - Ad ies (N, M) 



(5.2) 



and then adding the full NLO cross section, calculated using Eqs. (IAl|) - (jA5|) . which also 
includes the qg channels. 

VI. PHENOMENOLOGICAL RESULTS 

We are now ready to present some numerical results for our resummed rapidity-dependent 
cross sections and spin asymmetries for W production at RHIC. We are mainly aiming 
at investigating the quantitative effects of threshold resummation. In how far the spin 
asymmetry for this process can provide information on the polarized quark and anti-quark 



distributions has amply been discussed in the literature 



12, 15, 26, 401 and is not the 



focus of this work. We will therefore choose just one set of polarized parton distributions 



of the proton, namely the NLO "GRSV standard" set of j4l|. For the unpolarized parton 
distributions we take the NLO ones of |42| throughout. Unless stated otherwise, we choose 
the factorization and renormalization scales as \x = My/- 

Figures ^ (a) and (b) show the unpolarized and single-spin polarized cross sections for 
iy + -boson production in v^S* = 500 GeV collisions at RHIC, as functions of the rapidity 
of the boson. We show the results at various levels of perturbation theory. The lower line 
in Fig. (a) displays the unpolarized cross section at LO. The solid line above is the NLO 
result. The lines referred to as "resummed (l)-(3)" are for the matched NLL resummed cross 
section. As one can see, they all lie a few per cent higher than the NLO cross section. For the 
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FIG. 1: (a) Differential cross section for LF + -production in unpolarized pp collision at RHIC for 
V5 = 500 GeV as a function of the rapidity of the produced boson, (b) Same, but when one of the 
proton beams is longitudinally polarized, (c) The corresponding single-spin asymmetry Al- For 
all three plots, "resummed(l)" denotes the threshold-resummed cross section without including 
subleading ln(N)/N terms and without the cos^Mln-^ term, while for "resummed(2)" the 
Cosine term is included. "Resummed(3)" is with the ln(N)/N terms. The curve labelled "NNLO 
exp." gives the NNLO expansion of the "resummed(l)" result. We also display the LO and NLO 
cross sections. 



two dotted lines, denoted as "resummed (1)" and "resummed (2)", the subleading pieces 
oc \n(N)/N have been neglected. For "resummed (2)" the cos^Mln^) term discussed 
after Eq. (J3.9|) is kept in the Mellin moment integrand, whereas for "resummed (1)" it is 
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FIG. 2: Same as Fig. ^ but for W~ production. 

ignored. This evidently leads to a negligible difference. The curve for "resummed (3)" 
then shows the effect of also including the subleading terms oc hx(N)/N. This leads to a 
further small increase of the predicted cross section. Finally, the remaining line, "NNLO 
exp.", represents the two-loop (NNLO) expansion of the "resummed (1)" result. It shows 
that the NNLO terms generated by resummation are still significant, but that orders beyond 
NNLO have a negligible effect. We note that we have checked that the NLO expansion of the 
resummed cross section reproduces the exact NLO cross section very precisely. This provides 
confidence that the logarithmic terms that are subject to resummation indeed dominate the 
cross section, so that it is sensible to resum them. We use the same line coding for the 
polarized cross section in Fig. ^ (b). 

Overall, we find rather moderate resummation effects at \/S = 500 GeV. In particular, 
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FIG. 3: Same as Fig. □ but for y/S = 200 GeV. 



resummation does not affect the rapidity dependence of the cross section much, in the sense 
that the shapes of the curves from NLO to the NLL resummed case are very similar. This was 
to be anticipated from our analytical results in Sec. Hill where we found that the resummation 
of the rapidity-dependent cross section closely follows that of the rapidity-integrated one (see 
also 32]). Figure^(c) shows the corresponding single-spin asymmetry, obtained by dividing 
the results in (a) and (b). The higher-order effects cancel almost entirely, and resummation 
becomes unimportant. 

In Fig. El we display the corresponding results for W~ production. Very similar quan- 
titative features are found, except of course that the cross sections are smaller and, in the 
polarized case, of opposite sign because of the different combinations of parton distributions 
involved. 
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FIG. 4: Same as Fig. |2 but for y/S = 200 GeV. 



Threshold resummation should become more important when r increases. We thus expect 
larger perturbative corrections at RHIC's lower center-of-mass energy of y/S = 200 GeV. 
In Figs. El and |U we repeat our calculations in Figs. [T] and El at this energy. Indeed, the 
resummation effects are much more significant, amounting to about 25% at mid rapidity. 
Also, since we are closer to threshold now, subleading \n(N)/N terms are less important 
than before. Of course, the cross sections are much smaller at = 200 GeV than at 
500 GeV. In fact, luminosity will need to be at least around the design luminosity of 320 /pb 
in order to have sufficient statistics for a good measurement. It was shown recently 15] 
that at this luminosity a statistical accuracy of 5% (9%) for the unpolarized W + (W~) 
cross sections should be achievable, and that the current theoretical uncertainty related to 
the parton distributions is much larger, about 25%. It may therefore be quite possible to 
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FIG. 5: (a) Scale dependence of the NLO and the "resummed(l)" differential cross sections for 
W + production in unpolarized pp collisions at RHIC at \/~S = 500 GeV. We have varied the 
factorization/renormalization scale \x between and 2Mw- Note that we have multiplied the 

NLO cross section by 0.8 for better visibility. The uppermost curves in each case are for the scale 
Mw/2. (b) Same for the singly-polarized cross section, (c) Scale dependence of the single-spin 
asymmetry A^, calculated from the resummed cross sections. 

constrain better even the unpolarized parton distributions at large x from measurements of 
W ± production at RHIC at = 200 GeV. Clearly, it will be crucial then to take into 
account the perturbative corrections we show in Figs. El and EJ since these are of similar size 
as the current parton distribution uncertainties. 

So far, we have always chosen the factorization and renormalization scales as fi = M\y. 
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As we discussed in Sec. IIVI threshold resummation is expected to lead to a decrease of the 
scale dependence of the cross sections. We finally examine this in Fig. where we vary the 
scales to Myy/2 and 2M\y in the NLO and the NLL "resummed(l)" cross sections for W + 
production at y/S = 500 GeV. Note that we have multiplied the NLO cross section by 0.8 
for better visibility. As expected, the scale dependence is improved after resummation. At 
NLL, it becomes very small, and it cancels almost entirely in the spin asymmetry. 



VII. CONCLUSIONS AND OUTLOOK 

We have performed a study of perturbative higher-order corrections for W ± production 
in singly-polarized pp collisions at RHIC. This process will be used at RHIC to learn about 
the spin-dependent u and d distributions of the proton, and about the corresponding anti- 
quark distributions. In this work we have dealt with the resummation of potentially large 
"threshold" logarithms that arise when the incoming partons have just sufficient energy to 
produce the W boson. We have performed the resummation to next-to-leading logarithmic 
accuracy. We have considered the resummation for the rapidity dependence of the cross 



sections, for which we have used a method developed in [3l|. We find that the resummation 
effect on the unpolarized and single-longitudinally polarized cross sections is rather moderate 
at RHIC's higher energy y/S = 500 GeV, but more significant at y/S = 200 GeV where one 
is closer to threshold. The shapes of the W rapidity distributions are rather unaffected 
by resummation. We believe that our study will be important in future extractions of the 
polarized, and even the unpolarized, u, u, d, d distributions from forthcoming RHIC data. 

So far, we have only considered the W ± as the observed final-state particles. Ultimately, 
for the comparison to future data, it will be necessary to extend our study to a calculation 
of the cross section for a single charged lepton, since this is the signal accessible at RHIC. 
This will be devoted to a later study. We believe, however, that the results presented here 
already give an excellent picture of the resummation effects to be expected when the W 
decay is fully taken into account. 

We finally note that the technique adopted in this work for treating the rapidity- 
dependence in the resummation for Drell-Yan type processes should also be useful for ex- 
tending the recent study j^J of the Drell-Yan process in pp collisions at the GSI to the 
rapidity-dependent case. 
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APPENDIX A 



The NLO terms of the partonic coefficient functions in Eq. (j2.1j) read 19 
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